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We describe the Riemann-Hilbert (RH) approach to computing the long-time response of a Fermi 
gas to a time-dependent perturbation. The approach maps the problem onto a non-commuting RH 
problem. The method is non-perturbative, quite general and can be used to compute the Fermi gas 
response in driven (out of equilibrium) as well as equilibrium systems. It has the appealing feature 
of working directly with scattering amplitudes defined at the Fermi surface rather than with the 
bare Hamiltonian. We illustrate the power of the method by rederiving standard results for the 
core-hole and open-line Greens functions for the equilibrium Fermi edge singularity (FES) problem. 
We then show that the case of the non-separable potential can be solved non-perturbatively with 
no more effort than for the separable case. We compute the corresponding results for a biased (non- 
equilibrium) model tunneling device, similar to those used in single photon detectors, in which a 
photon absorption process can significantly change the conductance of the barrier. For times much 
larger than the inverse bias across the device, the response of the Fermi gases in the two electrodes 
shows that the equilibrium Fermi edge singularity is smoothed, shifted in frequency and becomes 
polarity-dependent. These results have a simple interpretation in terms of known results for the 
equilibrium case but with (in general complex-valued) combinations of elements of the scattering 
matrix replacing the equilibrium phase shifts. We also consider the shot noise spectrum of a tunnel 
junction subject to a time-dependent bias and demonstrate that the calculation is essentially the 
same as for the FES problem. For the case of a periodically driven device we show that the noise 
spectrum for the Coherent States of Alternating Current found in QJ can be easily obtained using 
this approach. 

PACS numbers: 72.10.Fk,73.23.Hk,73.40.Rw 



I. INTRODUCTION 

An approach to the study of the quantum statistics of 
an arbitrary single-particle observable in a Fermi gas 
has been recently described in We refer to it as the 
RH approach, as it reduces the calculation of a determi- 
nant describing the quantum statistics of an observable 
to the solution of an (in general) non-Abclian Riemann- 
Hilbert (RH) problem. The relation between such de- 
terminants and RH problems has been known for a long 
time, and has been used extensively in studies of quan- 
tum inverse scattering problems Q. As a result, a lot is 
known about the non-abelian RH problem Q , and much 
of this can be taken over directly to the study of the 
quantum statistics of Fermi gases. 

The RH approach gives an expression for the distribu- 
tion function of an observable in a Fermi gas perturbed 
by a time-dependent potential. To illustrate the method, 
one of us used it to prove a long-standing conjecture, first 
stated in || , that the two sources of shot noise in a biased 
point contact, namely fluctuations in the number of at- 
tempts to tunnel through the barrier and fluctuations in 
the number of reflections, are statistically independent 
0- We have also used the method to study how non- 
equilibrium effects alter the Fermi Edge Singularity in a 
tunnel junction [||. 

The response of a Fermi gas to a time-dependent per- 
turbing potential is a central problem in condensed mat- 
ter physics. It has been tackled in many different con- 
texts often with different approaches. For systems out of 



equilibrium, such as quantum pumps, perturbative ap- 
proaches, based on the Keldysh formalism, have been 
used, while for systems in equilibrium it has been possible 
to find exact solutions in some limiting cases by solving 
the equations of motion directly 0, 0, 0, E3 ■ O ne of the 
advantages of the RH method is that it applies equally 
to all such problems and therefore offers the prospect of 
a unified approach to computing the time-dependent re- 
sponse of all observables in Fermi gases. 

Setting up the description of a problem in the RH 
framework is quite straightforward. Given the solution of 
the single-particle scattering problem, the response of the 
Fermi gas reduces to the computation of a determinant 
of an operator taken over single-particle states occupied 
in the initial configuration. (The generalization of the 
method to the more general case in which the initial state 
is given in terms of a density matrix rather than a single 
quantum state should be possible but has not yet been 
formulated.) The evaluation of this determinant then re- 
duces to the solution of a Riemann Hilbert problem. The 
solution is in general a matrix-valued function analytic 
everywhere except across a cut, along which the function 
is discontinuous. The discontinuity is fixed by the driving 
force or perturbation acting upon the system Q. From 
the point of view of the Keldysh formalism the method 
performs a non-trivial re-summation of all relevant dia- 
grams with the help of the solution of the corresponding 
RH problem. In the abelian case, when the discontinu- 
ity function commutes with itself at all points along the 
cut, the solution is given in terms of an integral. The 
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classic solution of the FES problem 0, 0] is the simplest 
example of this solution. In the non-abelian case, the 
solution to the RH problem is not known in general, al- 
though asymptotic solutions exist. These are valid for 
response frequencies small compared to those present in 
the discontinuity function. 

Here wc explain the RH approach in some detail. To il- 
lustrate the power of the method we start by showing how 
the solution of the equilibrium FES problem [H, 0, is 
derived. We then show the generalization of this problem 
to include the case where the 'impurity' potential mixes 
scattering states of the unperturbed problem — the case 
of a non-separable potential — and deal explicitly with the 
case when the impurity potential gives rise to a bound 
state. This problem was treated initially in [ToL IT^ , in 
a calculation which solved directly the Dyson-like equa- 
tion for the appropriate Green's functions. In the RH 
formulation of this problem the discontinuity function, 
although matrix-valued, is constant and commutes with 
itself. As a consequence, the solution to the RH prob- 
lem is trivial to derive and yields the standard results 
of 0, ^3 with no more work than for the case sep- 
arable potential case. We show how these results are 
changed in a non-equilibrium situation. In both the equi- 
librium and non-equilibrium cases, we compute both the 
core-hole Green's function reported in Q and the open 
line contribution. Finally we show how the states which 
minimize the shot noise in a periodically driven quan- 
tum pump — the so-called Coherent States of Alternating 
Current (CSAC's) Q 

— can be described using the RH 

method. 



II. PERTURBING THE FERMI GAS 

We consider a system in which particles impinge upon 
a localized potential (see FigQJ. The potential is time- 
independent for all times t < to and t > tf. For times 
to < t < tf the potential varies. We take our basis to be 
the eigen states of the system with the localised potential 
at its value. The states are labeled by their single-particle 
energy e (fr = 1) and a channel index i. We will consider 
the corresponding annihilation operator, a^, as the i'th 
component of the vector a e . The Hamiltonian of the 
system is then 

Ho = ]Tea+a £ . (1) 

e 

Here M(i, e, e') = when t < or t > tf. (In the follow- 
ing,^ for any operator O, we will denote by O the matrix 
of O taken between the single-particle basis states.) 

We will be interested in the total effect of the pertur- 
bation, ie what is the final state of the system for t > tf 
given the initial state at t = 0. This requires a knowl- 
edge of the effect on the initial many-body state of the 




FIG. 1: Schematic illustration of the generic scattering prob- 
lem. Particles impinge on a localized potential which is in- 
dependent of time for times t < to and t > tf. For times 
to < t < tf the potential varies as a function of time. The 
instantaneous value of the potential at time t is shown as the 
dashed curve and the difference measured with respect to the 
time-independent potential is denoted by M(t). When the 
potential varies slowly (condition |KJ , the effect of the scat- 
tering potential on an incident partial wave (shown schemat- 
ically as (j)) is simply multiplication by the scattering matrix 
corresponding to the instantaneous value of M(t). 

time-development operator U(tf), where 

i^=H(t)U(t), U(0) = 1. (2) 

Because the Hamiltonian H(t) in is quadratic, the 
effect of U{tf) is fully characterized by its effect on the 
set of single-particle scattering states, af e , |): 

U(t J )k+\} = J2e- i ^a(e,e')kt\), (3) 

€ 

where cr(e,e') is some unitary N x N and |) is the true 
vacuum with no particles in the system [1^ . 

When computing the response of the Fermi gas to the 
time-dependent potential, we will need to compute ex- 
pectation values of the type: 

xr = (o\R\o)- (4) 

Here |0) is the state of the Fermi gas before the perturba- 
tion is applied and R is an operator (or operator product) 
related to an obvcrvable of interest. In general, the R 
in will involve the time-development operator U(tf). 
For example, in the case of the shot noise spectrum of 
a tunneling barrier 0, the interest is in the statistics of 
the charged transferred from one electrode to the other. 
If Qi is the charge in the first electrode then the expec- 
tation value of R = U\t f )e iX ^ U^^e' 1 ^ 1 , yields the 
generating function for moments of the distribution of 
charge transferred out of channel 1 (into channel 2) dur- 
ing the period between t = and t — tf. In the case of 
the FES problem @, the core-hole Green's function (see 
below) is related to the overlap (0\Uo(tf)U(tf)\0) where 
Uo(tf) is the time-development operator for Hq in 
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This overlap is an expectation value of the type J3J with 
R = U^t f )U(t f ). 

The effect of U{tf) acting on the single-particle states 
of the basis is given by the unitary matrix a(e, e') defined 
in The matrix er(e, e') can be related to the scattering 
matrix, S(t, E), for a particle with energy E evaluated on 
the instantaneous value of the potential M(t) in JQ). This 
reflects the fact that S(t, E) encodes all the information 
in the potential variations M(t). This relation will be 
complicated in general. However, when 



^dS- 1 dS 
h —dE <<1 ' 
the relation between a and S becomes simple: 



~>uv 



1 



dtS lv {t,E)e^-^\ 



(5) 



(G) 



where E = (e + e')/2 and Vi is the density of states in 
channel i. The result © shows that the total scatter- 
ing amplitude from state k in channel n to k' in channel 
n' is just the Fourier transform of the scattering matrix 
S(t) evaluated on the instantaneous value of the poten- 
tial. This result is well-known. It was used implicitly 
to solve the FES problem in the presence of a separa- 
ble potential ^lj and is often used to simplify studies of 
quantum pumps. 

A brief derivation of the condition 10 is given in Q . 
The condition can be understood heuristicaUy as follows 
(see also pj). We consider the incoming wave-packet to 
be a partial wave in channel n of the basis, in which S 
is diagonal before the perturbation is switched on. After 
impinging on the potential, the partial waves scattered 
from channel n into channel n' will take a time of order 
of the corresponding Wigner delay time to pass out of 
the region where the potential acts. The condition (J3J) 
is equivalent to the requirement that the scattering ma- 
trix does not change significantly during this delay time. 
If this condition is satisfied, the relation © also has a 
simple interpretation. A wave with energy e' = k'v in 
channel n' incident at time t on the potential (which is 
assumed to be localized around the origin) will have am- 
plitude at the origin proportional to e~ lk vt where v is the 
velocity of the incident wave. Waves will emanate from 
the source at the origin with amplitude S nn '(t)e~ lk vt in 
channel n. If dispersion effects are small this will lead to 
a waveform S nn '(t~(r/v'))e~' lk v ( t -( x / v )). Decomposing 
this into waves of the form e lk ( x ~ v t) as i _> ^ gives the 
result ©. The normalization factor 1/ \Jviv[ is included 
so that in the case where the incoming flux (~ v\ip n \ ) is 
totally scattered into channel n' the scattering amplitude 
is 1. 



A. Fermi Sea at T = 

The calculation of the response of the Fermi gas to 
the time-dependent perturbation reduces to the compu- 



tation of the expectation value XR = (0|-R|0) in In 
the following we will assume that |0) is a single Slater de- 
terminant. In this case the computation of XR requires 
the evaluation of a single determinant : 



(0|#|0) = det' | R |, 
where the elements of R are given by: 

Rivie,^) = (|6 ie i?4 e ,|). 



(7) 



(8) 



and where the prime on the det indicates that the deter- 
minant is taken only over states occupied in |0). 

When the initial Slater determinant |0) corresponds 
to a filled Fermi sea, it is useful to introduce the Fermi 
distribution as an operator with elements: 



8e, e '6ur6(-(e - fx)). 



(9) 



Here fj, is the chemical potential, which we take to be 
zero. (For the non-equilibrium problems discussed later 
the chemical potential can vary according to the chan- 
nel index. However, a time-dependent change of basis 
replaces the differing chemical potentials in the different 
electrodes by an additional time-dependence in the scat- 
tering matrix S(t), so that there is no loss of generality 
in assuming \x = 0.) In a block notation that separates 
the states with positive and negative energies / and R 
become: 



/ = 



1 




R = 



Rn R\2 

i?21 R22 



and 



i-f+fR 



Rn Rv. 







1 



It then follows that 



XR 



det' I R n \= det|l - / + fR\ 



(10) 



(11) 



(12) 



where det is now a determinant taken over all states in 
the basis. 

Expressing xr as the determinant of 1 — / + f R taken 
over all states in the basis, allows us to write 

log xr = Tr (/ In R) +Tr (In (1 - / + fR) - f In R) (13) 

In i|13fl . we have added and subtracted the term 
l°gXfl' = Tr(/lni?). This term consists of the diago- 
nal elements of In R summed over all occupied states in 
|0) and gives the contributions linear in tf. It often has 
a simple physical interpretation. In the FES problem it 
yields the threshold shift (or change in the ground state 
energy of the Fermi gas after the core hole is created), 
while in the shot noise spectrum of the tunneling barrier 
it can be shown to be related to the average transfer of 
charge across the barrier (the Brouwer formula 0, El ) • 
The second term in <|13|) . 



log X ^ =Tr(ln(l 



f + fR)-flnR). 



(14) 
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accounts for all the non-trivial effects associated with 
excitations close to the Fermi surface induced by the 
perturbation. (States far from the Fermi energy, when 
/ = 1 or / = 0, make no contribution to this term. As 
a result there are no problems associated with effects of 
the band edge or short time cutoff when computing this 
term.) In the case of the FES problem it describes the 
line-shape, while in the shot noise spectrum it gives all 
the higher moments of the charge transfer distribution. 
When computing this term, we will later switch to the 
time-representation in which 



Rlel'e' 



1 



dtR u >(t,E)e 



i(e— e )t 



(15) 



with (as in|SJ E = (e + e')/2. R will normally involve 
a or some simple combination of a with itself and its 
inverse. Provided the condition 10 is satisfied, we will be 
able to evaluate a by ignoring the dependence of S(t, E) 
on E. (This E dependence is not important as states 

(2) 

far from the Fermi energy do not contribute to logx^ ■) 

(2) 

As a result the term logX/? will depend only on the 
time-dependence of the scattering matrix evaluated at 
the Fermi energy, which we will denote by S(t). 



B. The Riemann-Hilbert Problem 

Computing the second term in l|13|) , log Xr ; is the cen- 
tral task in the evaluation of the response of the Fermi 
gas. The non-trivial part of this is finding the inverse of 
(1 — / + f R) which can then be used in an integral repre- 
sentation for its logarithm. This inverse can be written 
in terms of the solution of an TV x TV matrix Riemann- 
Hilbert problem, where TV is the length of the vector a e , 
ie the number of channels in the problem (sec0. 

A standard procedure for representing the logarithm 
of an infinite matrix, such as the one on the right hand 
side of lfH|) . introduces a A-dependence for R via ^5| 



R{\) = exp(Alogi?) 



(16) 



and then uses an integral over A to represent the loga- 
rithm: 



(2) 



dATr 



((i-Z + ji?)- 1 /-/^ 1 )^ 



(17) 
in 



(The A-dependencc of R introduced in i|16|) is assumed 
(|17fl although not written explicitly.) 

To compute the trace in (|17|) , we switch to a time rep 
resentation in which a quantity A becomes: 



A n ,{t,t') = - / vide / u v de'A lv (e,e') (18) 

Now, the Fermi distribution @ is no longer diagonal: 



2-Kt-t' -MO' 



(19) 



However, as we can neglect the E dependence of 
R(X,t 7 E) and S(t,E) (see discussion after ^J, R and S 
are now diagonal in t and equal to R(X, t, 0)S(t — t') and 
S(t, 0)5(t — t 1 ) respectively. In the time-representation, 
the product of two quantities requires matrix multiplica- 
tion in the space of scattering channels together with an 
integral over the intermediate time coordinate. Where 
one of the quantities in the product is diagonal in t (for 
example S), the integral over the intermediate time co- 
ordinate is of course trivial and the product reduces to 
the simple matrix multiplication in the channel space. 
The Tr now becomes a trace over the scattering chan- 
nels, which we denote by tr, and an integral over the 
time coordinate, so that 



(2) _ 



d\ / dttr 



((l-f + fR)-\f-fR-i)^. 



(20) 

Here, when computing the diagonal (equal time) ele- 
ments of O, one should take lim t ^ t / 0(t, t'). If A and 
B are diagonal in the time representation, it follows that 



dttr[A, f]B 



dt lim tr — 

t'-»t 2tt 



A(t)-A(t<) 
t - V + iO 



B{t') 
(21) 



which is a result we use later. 

The quantity (1 — / + JR)^ 1 can be written in terms 
of the function Y(t), which is a matrix in the chan- 
nel space and which solves an auxiliary Riemann-Hilbert 
problem. Y(t) should be analytic everywhere in the com- 
plex t— plane except on the interval [0, t/] on the real axis 
along which it satisfies 

y_(t)y + " 1 {t) = R(t) where Y± = Y(t±i0). (22) 



In addition Y should satisfy 



Y -> 1 



when 



1*1 



(23) 



These analytic properties together with (|19fl yield the 
useful identities 



fY-f 
fY+f 



fY- 

y+f 



(24) 



Using these relations (and assuming that F _1 is also an- 
alytic everywhere except along the cut), it is then easy 
to verify that 



(1 - / + fR)- 1 = Y + ((1 - f)Y+ 1 + fYl 



(25) 



As an aside, we note that (1 — / + fR) -1 is the solution 
to a singular integral equation, with / playing the role of 
the singular kernel of the Cauchy type. It is well known 
that such integral equations can be solved using Carle- 
man's method, which writes the solution in terms of an 
analytic function satisfying a Riemann-Hilbert problem 
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[16J . In the one-channel the corresponding singular inte- 
gral equation for the case when R(t) is constant between 
t = and t = tf is the problem solved in (3, 0] when 
describing the equilibrium FES. 

Inserting (|25|l into (|20|l and using (|21|l we obtain 

The integral over t is over all times. However, as dR/dX 
normally vanishes for t > tf and t < 0, one often only 
needs to integrate from to tf. 

Equations l|2t)|) and (|13|) map the characterisation of 
the response \R m Q onto an integral involving the so- 
lution, Y(t), of a Riemann-Hilbert (RH) problem 1)221 and 
I2SJ- An appealing feature of this formulation is that these 
formulas apply for any choice of variable \R provided 
that the Fermi gas (or gases) is initially in its ground 
state and apply for many non-equilibrium cases as well. 
If R(t) commutes with itself at different values of t, the 
solution of the RH problem can be written in closed form. 
Although there is no solution for the general case, a lot 
is known about such non-commuting problems including 
some asymptotic solutions valid when tj 1 is much smaller 
than any characteristic frequency in R(t) 0. 

III. FERMI EDGE SINGULARITY 

In this section we show how all the known results for 
the equilibrium Fermi Edge Singularity (FES) follow di- 
rectly from the formula i|26|) . Within our formalism the 
case of non-separable channels considered in [Tol | and 
again in is n0 more complicated than the separable 
case. We will then discuss how these results are changed 
in the nonequilibrium case. 



A. Equilibrium 

The FES problem was first considered in the context 
of the X-ray absorption spectrum of a metal Q ■ When 
a photon creates a core hole in a metal, the Fermi gas is 
affected by the potential of the core hole leading to the 
excitation of particle-hole pairs. The absorption line ex- 
pected in the absence of the Fermi gas becomes a thresh- 
old with a singularity in the absorption spectrum as a 
function of lo — loq > 0: 

/(w)~|w-w |- a , (27) 

where ojq is the threshold frequency for absorption. It 
turns out that similar singularities are seen in the distri- 
bution of energy absorbed by the Fermi gas in response 
to any rapid change in potential and not just in x-ray 
absorption experiments. For example, the consequences 
of the FES are also seen in a tunnel junction. As the 
energy absorbed by the Fermi gas when switching, is an 



important characteristic of the device, establishing how 
the FES changes in such tunneling devices is important 
for understanding fluctuations in energy transfer across 
such devices. The FES is also thought to be related to 
the apparent absence of detailed balance in Random Tele- 
graph Signals [l7j . 

The Hamiltonian for the photon absorption experiment 
is H 

H = e<$b + J2 £ a+a e + k\V (e, t')k t , ffi + H x (28) 

e e'e 

with the operators a as in P and V(e, e') an TV x N 
matrix. The operator w is the creation operator corre- 
sponding to the core state and the coupling to the X-ray 
field is described semiclassically by 

Hx = '« ' k \be lut + h.c. 

e 

= Xe iuJt (29) 

The absorption spectrum is proportional to the real part 
of the Fourier transform of the response function 

S(tf) = (0\T{X(t f )X(0)}\0), (30) 

with T the time-ordering operator. S(tf) can be com- 
puted from the core-hole Green's function 

G(t f ) = (0\T{bHt f )b(0)}\0) (31) 
and the function 

F(t f ) = ^(0|T{St(i / )(W e Va e ,(i/))}(W e -at(0))S(0)|0). 

. (32) 

Coventionally a minus sign is included in the definition of 
F and G. However, as wc will only deal with the absorp- 
tion case here and take tf > 0, it is easier to work from 
these definition. We have also left out the conventional 
factor of i in the definitions of these Green's functions as 
in 0. 

The calculation of F and G reduces to a one-body 
scattering problem 0, [^. As far as the Fermi gas is 
concerned the role of the core hole is to switch on the 
scattering potential V(e, e') at time and switch it off 
again at tf. As such, the problem is clearly in the form 
of m with M(i,e,e') = for t > tf and t < and 
M(t,e,e') = V(e,e') for < t < tf. The corresponding 
scattering matrix S(t, E) switches between the identity 
when the core-hole is absent and some constant value 
S e (E) when it is present. The asymptotic behavior of 
the response at large tf (when lo — lo <C ^q 1 ) is deter- 
mined by states with energies close to the Fermi surface. 
For these states we assume that the variation of S e (E) 
with E can be neglected so that the condition JSJ) is sat- 
isfied. (The limit £o*/ 3> 1 is the one considered in 0.) 

The calculation of G is one of the simplest calculations 
within the RH approach. G is the expectation value of 
the operator R in Q with: 

R=ul(tf)U(tf). (33) 
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As the matrix elements of U(tf) are just e ietf a ee /, it 
follows from that [HI 

G(t f ) = e leot 'det' | a | (34) 

while from (pTI 1151 and 116}! 

R(t) = S{t) and R(X,t) = expAlog5(i). (35) 

The RH problem ^undfEft reduces to 

Y-(t)Y+\t) = S x {t) Y -> 1 when |i| -> oo. (36) 

When the matrix S is constant between and tf, we 
will denote its value by S e . In the single-channel case, 
S e = e 2tS and the RH problem is solved by [l^ 



Y{z) = cxp 



1 

2tti 



In 



AlogS £ 



(37) 



(This solution was used implicitly in the original solution 
to the single channel problem of @-) In fact the result 
(|37|l solves the RH problem even where the problem is not 
separable provided that the matrices S(t) evaluated at 
different times t with < t < tf commute. (This can be 
checked by direct substitution into H36|> .) We insert Y(z) 
and R(X,t) into (|26[l . The integral over t runs between 
and tf where logi?(A) is non-zero. Inserting into l|13|) 
and including the factor of e le ° tf yields 



logXi? 



ie' t 



logi£o*/ 



(38) 



where e' = eo + J2 e <o <H e )/ 7r ^ / ( e )j with i/(e) the density 
of states, is the shifted energy of the core-hole in the 
presence of the Fermi gas. (The form for the difference 
between eo and e is usually attributed to Fumi |Tii.l2flj.) 
Close to the branch points of Y at and tf, we cut the 
integrals off at i^Q 1 and tf + i^Q 1 where £o is an energy 
of order the band width. Eq l|38|) gives the well-known 
result for the long-time asymptotic behavior of G Q: 



G(t f ) ~ (itotr^f, 



a 



(39) 



To compute the function F(tf) in 1|32|) is slightly more 
involved, although the underlying RH problem is the 
same. As already mentioned, the role of the core hole 
is to switch on the potential at t = and switch it off 
again at tf, so that F(tf) can be written (writing out the 
channel indices explicitly) 

F(tf) = E WZ(0\f(ie,M)\0)W Ve , (40) 
f(ie,i'e') = Ut(t f )a ie U(tf)4e>- («) 



In the basis of the scattering states aj, a ,\) the matrix 
elements of this operator are easily shown to be given in 
terms of cr in by 



r(ie, i'e')j 



jaj a' ^ let' e' 



(42) 



Using J7J) and gUJ) we find 

F(t f ) = e itot fdet' \ Ca - \h){g\ | 
Here C = C(tf) is the number: 



ie^iei'e'Wi'e' 



(43) 



(44) 



while 



\h) 



{g\ = 



ja \i'e' / 



The expression l|43|) is now in the form |J7J|. We could 
attempt to solve the corresponding RH problem l]_!2l) 
and (|23|l as before, although the relation between the 
corresponding operator R{t) and S is no longer simple. 
However, it is easier to simplify l|46|) by factoring out 
G(t f )=e ieot fdet'\a\: 



with 



F(tf) = CG(t f ) det | 1 - C~ 1 0\h){g\ I (46) 



O = (1 - / + fa)- 1 !. (47) 



We have used Q12[l to put l|46() in the form of the deter- 
minant over all states in the basis. Using the identity 
det|l - C- 1 0\h){g\\ = 1 - C- 1 {g\0\h) we obtain 



F(t f ) = G(tf)(C-(g\0\h)). 



(48) 



As C — * for large tf (with a functional form which 
depends on assumptions about the density of states at 
the band edge) , the response is determined by the second 
term. 

The function Y(z) with A = 1 in (|5B1 and l23|l can be 

used to compute F(tf). In the time-representation 



(l-f + fa)- 1 f = Y + fY- 



(49) 



Writing F(tf) = L(tf)G(tf) (L is usually referred to as 
the open-line contribution), we find 



L = -E / dede ' \ dt 1 dt 2 W* e e ieit:L - t f' ) x 



v v e " ie ' t2 W Vel (50) 



Taking Wi to independent of e (we arc assuming that the 
long tf behavior is determined by states with energies 
within ~ 1/iy of the Fermi surface), this simplifies to 
give 



1 



(51) 
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The functions Y_ and Y+ are evaluated at t = and i = 
tj . This prescription is equivalent to the imaginary time 
cutoff used to derive and used in 0, . Strictly, 
the discontinuities in S at t = and t = tf should be 
thought of as the limit of a fast switching process, in 
which S starts to change at t = and reaches its new 
value S e after a short time. Similarly at t — tf, S starts 
to change back from S e to its unperturbed value. (The 
corrections associated with a more realistic model of a 
non-instantaneous switching process were considered for 
a related problem in [2l]].) In the single channel case 
we can insert the explicit form for Y given by (|37|l . and 
recover the standard results 



1 



1 



1 



1 



it f (i£ t f y s M 2 - 2S /* 



(52) 



When the potential M(t) in (Q is strong enough for 
a bound state of the Fermi gas electrons to form below 
the bottom of the band, the results for G(tf) and F(tf) 
are no longer correct. Provided the perturbing poten- 
tial is constant between and t / , the effect of the bound 
state can be taken into account explicitly as explained 
in the appendix. The result for G(t/) given by HA2jl. 
and for F(tf) given by l|A12jl . have two main contribu- 
tions. After taking the Fourier transform to obtain the 
absorption spectrum, the first corresponds to having the 
bound state occupied and leads to the absolute thresh- 
old for absorption. The second term relates to scattering 
processes in which the bound state is always empty and 
leads to a subsidiary threshold at Ei, above the first in 
the absorption spectrum. 

The results (j^J), ljB"2*|l. l|A13j) are of course very well- 
known However, none of the key formulas (|37l 1511 
and IA12T) require that the scattering matrix S should be 
diagonal in the channel indices. Provided that S(t) com- 
mutes with itself at different times, the results are valid 
for arbitrary channel number. We can therefore use the 
function Y{z) given by l|37|) to compute the correspond- 
ing results for the case of a non-separable potential just 
as easily as in the separable case. In the absence of bound 
states, one obtains with e' = + S e <o ^c( e )/ 7rI/ c( e ) : 



G(t f ) = exp(ie' t f ) (i£ t f 



Htf) 



t f 



(53) 



with defined as the phase shift modulo n in channel £ 
on the interval [—ir/2, ir/2] (see discussion after IA7|) . We 
then obtain the gene ralizations to the non-separable case 
of the results of [111 for G(t f ) and F(t f ). We find 



G{t f )=G{t s ){l + A B ) 



(54) 



where G(tf) is the contribution of the scattering states 
given by the expression for G in (|53|l with phase shifts 
given by 5^, while 



A B ~ e- lE ^ J2 k| 2 e _2i5< cxp ( ~5 C lnifot/ 

C 



(55) 

Here rj^ = J^i V^i u * f(C)i an d the u,; are the bound state 
wavefunction coefficients given in l|Alfl. In the presence 
of the bound state, the function F(tf) ~ Fo(tf) + Fb(tf) 
with 



Fo(t f ) 
F b (t f ) 



G(tf)L(t f ) 

e - iEBt fG(t f )\u- W| 2 , 



(56) 



where L(tf) is the scattering state contribution to L(tf) 
given by the expression in ij53*|) , using the phase shifts 5^ . 

The formulas (|56|l and l|53|l are the natural generaliza- 
tions of the single channel result and have exactly the 
same interpretation as was given originally in l22j| . 
We repeat this briefly here as the results for the non- 
equilibrium case (given in the next section) can also be 
understood heuristically on a similar basis but with the 
phase shifts becoming complex. The exponents {S^/tt) 2 
and (<5^/7r±l) 2 are, according to the Friedel sum rule, the 
square of the net charge that needs to move in to or away 
from the origin in order to screen the core hole potential. 
For G(tf) this is S^/n, while for F(tf) it is (5^ — tt)/tt if 
the photoclcctron inserted at the origin is in the £ chan- 
nel and otherwise. If there is an occupied bound 
state after absorption of the photon, the respective val- 
ues become (5q+tt)/tt and 5q/-k, as now the Fermi gas has 
to provide the additional electron which ends up in the 
bound state. The form t~ n is just the decay with time 
of the overlap of the wavefunction of the Fermi gas at 
t = tf and the one describing the system created at t = 
in which (with respect to the ground state in the pres- 
ence of the core hole) there is an excess charge n = —5/tt 
at the origin. That it vanishes as t — > oo, is the orthogo- 
nality catastrophe described by Anderson [23T |. 



which are the results obtained perturbatively in [l2^ . 
Here the eigen values of the matrix S e (see after I36f) 
are written as e 2lSc - . S e has eigen vectors /(C)i and 

The perturbing potential, characterised by scattering 
matrix S e , can be strong enough to lead to a bound state 
with wavefunction given by (|A4|) . In the presence of a 

bound state(s) we take the eigenvalues of S e to be e l2S( 



B. Non-equilibrium Effects 

The experimental and technological interest in the out- 
of-equilibrium response of coupled Fermi systems has 
grown as electronic devices have shrunk. Examples in- 
clude structured quantum dots, like the single electron 
transistor or the single photon detector |24|, and quan- 
tum point contacts. The non-equilibrium Fermi Edge 
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FIG. 2: Energy levels in an idealized device to demonstrate 
the out-of-equilibrium FES. The scattering potential for elec- 
trons is characterized via the 2x2 matrix, S(e), connecting 
scattering states in the two wires for particles with energy e. 
S — S 9 or S e depending on whether the defect is in its ground 
(g) or excited (e) state (with excitation energy _Eo). S 9 is the 
identity matrix and S e is an arbitrary unitary matrix. Sfi 
and Sf 2 correspond to the reflection and transmission ampli- 
tudes respectively. We will refer to the device operating as 
illustrated here, with a negative potential —V (V > 0) applied 
to the left electrode, as the forward- biased case. 



Singularity (nFES) will characterize the energy absorbed 
by the coupled Fermi gases in a rapid switching process 
in such devices. The nFES should help explain, for ex- 
ample, measurements of random telegraph signals (RTS). 
In these experiments, a two-level system (TLS) couples 
to the source-drain current flowing in the channel of a 
MOSFET (the TLS resides in the insulating oxide layer 
The RTS relates to the 'random' switching of the 
TLS between its ground and excited states. The ratio 
between the times the TLS spends in the excited and 
ground states is measured experimentally. In equilibrium 
this ratio is fixed by detailed balance, and the deviations 
from this have been attributed to non-equilibrium effects 

m. 

As one of the simplest non-trivial many-body effects, 
the FES is also a natural point to start, when looking for 
a description of non-equilibrium effects in many-electron 
systems. Perhaps surprisingly, given its conceptual sim- 
plicity, the nFES has not attracted as much attention as 
more difficult non-equilibrium problems like the Kondo 
effect, to which it is known to be related. (The Kondo 
effect can be thought of as a sequence of FES's associated 
with each flipping of the localized moment H^|.) 

In we reported results for G(tf) for a two-channel 
problem, which modelled a system with two electrodes 
separated by a barrier. The transmission of the bar- 
rier depends on the state of a two-level system inside 
the barrier, see Figure [21 with the transition between 
the two-levels assumed to be dipolar. The real part of 
the Fourier transform of the function G(tf) gives the ab- 
sorption spectrum for the device. The non-equilibrium 
effects predicted in should be visible in the voltage de- 



pendence of the absorption line-shape of devices like the 
single-photon detector of [24[ and [23 . In [24[ a quantum 
dot in the quantum Hall regime is coupled via tunneling 
barriers to two electrodes on either side of the dot. For 
magnetic fields in the range 3.4-4.2T, the conductance 
through the dot can change from zero to around 0.3 e 2 /h 
when a photon is absorbed via cyclotron resonance in 
the dot. From the perspective of the two electrodes, the 
dot behaves as a tunneling barrier which allows tunneling 
only in its excited state. The absorption of the photon 
and the subsequent separation of the hole (which moves 
into the v = 1 ring on the outer part of the dot) and the 
particle (which 'falls' into center of the dot at v = 2) is 
rapid, while the response of the conduction electrons in 
the two electrodes is slow and will show effects character- 
istic of the FES. In the device of |27|. a electron trapped 
in a dot underneath an electron channel gives rise to a 
potential which closes off a conducting channel. When a 
photon is absorbed, the photo-exicted hole can recombinc 
with the electron in the trap, the potential of the electron 
disappears and the channel opens. Again the conduction 
electrons on the two sides of the channel, 'see' the sud- 
den reduction of a tunneling barrier on absorption of a 
photon. 

The main result reported in [(| was that the ND for- 
mula describing the form of the FES and threshold shift 
(Fumi's theorem [THH) generalized in a simple way to 
the non-equilibrium case. For time scales i < tj, the 
phase shifts which appear in G(tf) are real and are given 
by the log of the eigen- values of the scattering matrix S e 
\wL Il2| ■ This simply reflects the fact that on these short 
time scales the response of the Fermi gas involves ex- 
citations with energies e > 7 which do not sense the 
non-equilibrium distribution function. On time scales 
t ^> tf, the equilibrium phase shifts in the two channels 
are replaced by 'complex' phase shifts given by log 5^ 
and log (1/5*22)* • The real part of these phase shifts de- 
scribes the scattering within each electrode, while the 
imaginary part describes the effect of scattering processes 
in which particles cross the barrier. One effect of the 
non-equilibrium operation of the device is to make the 
scattering between the different electrodes effectively in- 
coherent. Here, we find that this interpretation extends 
also for the function F(tf). 

We show the key steps in the derivation of G(tf), em- 
phasizing the relationship with the equilibrium results, 
and report the results for F(tf) including the role of pos- 
sible bound states. Since the initial state involves a filled 
Fermi sea in both channels (left and right electrodes), 
the RH formulation of this non-equilibrium problem is 
the same as for the equilibrium case. The bias across the 
tunnel junction means only that the chemical potentials 
arc different in the two electrodes. One way of handling 
this difference is to introduce a gauge transformation act- 
ing only on the basis states in the left electrode: 

a(e) -> a(e,t) =e +iL ti v ^ dT a(e) (57) 
S(t) -> S(t) = e +lL -fo v ^ dT S(t)e~ lL -fo v ^ dT (58) 
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where L is the diagonal matrix with L\\ = 1 and 
L22 = 0. The effect of this transformation on states 
in the left electrode is to set e — > e — V(t), so that the 
chemical potential in the left electrode becomes equal 
to that in the right electrode (taken to be zero as be- 
fore). For the constant bias case, the transformation 
gives a(t) -> a(t) = e+ iLV 'a. 

(2) 

The functions logx^j and L(tf) for the nFES case are 
still given by l|26|) and (|50[1 . However, the RH problem 
satisfied by the function Y(z) is different: In 1|36|) S* A (t) 
picks up an additional time-dependence from the gauge 
transformation (|58f) . which leads to two important dif- 
ferences to the equilibrium case. Firstly, the function 
e lVt introduces a new characteristic energy scale, V. If 
the function S(t) has Fourier components with freqen- 
cies u> ^> V, the response will be dominated by states 
with energies |e| ^ V and will be insensitive to the 
non-equilibrium nature of the distribution which only be- 
comes apparent on the energy scale V. If S(t) only varies 
at frequencies u> <C V the response will come from states 
with energies e <C V and will normally be significantly 
different from what happens in equilibrium. 

The second main difference following from the addi- 
tional time-dependence of S(t) relates to the case when 
between t = and t — tf the scattering matrix (be- 
fore the gauge transformation) is constant and equal to 
S e . In this case it is now no longer possible to solve 
the RH problem with a function of the form (|37[l . Al- 
though this form satisfies formally the jump condition, 
Y—(t)Yr(t) = S*(t), the corresponding function Y(z) 
is not well-defined for large z if S e is not diagonal. The 
off-diagonal elements of S e contain factors e ±lVt so that 
there is an essential singularity at z — > 00 in Y{z) de- 
fined by H37|) . an d it no longer satisfies the condition 
Y — ► 1. This problem is clearly apparent in the RH 
formulation we have presented. It was much less clear 
in previous attempts to extend the ND method to the 
non-equilibrium case and may explain why these failed 
p8j . It is also interesting to note that for tf <C 1/V, wc 
can expand the function e tVt up to linear order in Vtf. 
Then S e = S e (V = 0) + CVt, there is no singularity at 
infinity for Y, and the form l|37() still works. 

In general there is no exact solution to the non- 
commuting RH problem Q. However, in the case rele- 
vant to the device shown in FigEl S(t) = e lLVt S e e~ iLVt 
for < t < tf, with S e constant, we can find an asymp- 
totically correct solution for the limit tf 3> 1/V relevant 
to the nFES 0, El • We will only consider the case where 
there is one channel in each electrode. As in the equilib- 
rium case (cf. l33|) 



i? = cxp (AlogS(t)). 



(59) 



In this case the solution for Y(t), valid for t ^> V 1 , is 
given for t < or t > tf by 



while above and below the cut, [0, £/], 

'1 -i{t A) 



Y+(t,X) 
Y-(t,\) 



I 

1 

+»?(*, A) 1 



lM*,A). (61) 



Here j(t, A) = Ru/Rn and r?(i,A) = R 2 i/Rn. The 
functions ip±(t, A) = tp(t ± iO, A), where tp(z, A) is given 
by 



l\) = cxp 



(xiT +X2T 3 )l0g- 



(62) 



with 



logRn/R* 22 , m log(R n Rh) 
xi(A) = — ; and X2(A) = 



47TI 



4iri 



(63) 

Here T3 is the third Pauli spin matrix and To is the iden- 
tity matrix. The derivation of (|61|1 follows that given 
in 0- The idea, which was explained in detail in the 
context of inverse scattering problems in , is to solve 
for a function W-^z) which satisfies the same jump con- 
dition as Y(z) but in a complex plane with additional 
cuts. For this problem, the additional cuts are parallel 
to the imaginary axis and run from the branch points at 
3 = and z = tf to infinity. The discontinuities in W(z) 
across the vertical cuts scale as e - '^ 2 '. If Y is approxi- 
mated by W, the errors in \ogxn defined in l|13|) are only 
0{1/Vtf), and can, in principle, be computed order by 
order in powers of (Vtf)^ 1 . 

The form for log xr for tf 3> V^ 1 is found by inserting 
(|61fl into (|26|) and (|13|l and computing the integrals over 
t and A as in the equilibrium case 0]: 

log X (tf,V) = -t(E -A(V))tf-f3'\og(tVt f )+D, (64) 

where A(y) is given by the non-equilibrium generaliza- 
tion of Fumi's theorem 0, [2(1 



A(V) = f 



2-7TZ ./o 







trlog (S*(E))_ dE r v log (Sf^E)) d£ 



(65) 
(66) 



The constant (3' is given by (cf. I53fl : 

•log (sfo \ 2 , (\og(i/ s% 2 y 



2-ni 



2ni 



The constant term D can be estimated by requiring that 
the form for logx l|65|l matches the equilibrium one at 
tf = V^ 1 (GEll valid for t f < V -1 . This constant gives 
the contribution from excitations with frequencies be- 
tween V and £o- This gives: 



D = [3log£ /V. 



(67) 



Y(t,X) = V(i,A), 



(60) 



The result for G(tf) can be seen as an adaptation of 
the equilibrium result. The real phase shifts (given by 
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FIG. 3: Spectral function Kexd^) computed from <B10t 
with ui in units of the bias voltage V. The spectra depend 
on Sf- L = ^/Re i2ai and Sh = VRe i2a2 , where R is the reflec- 
tion probability. The curve marked + (— ) refers to the case in 
which electrode 1 (2) is at the higher chemical potential. Also 
shown is the corresponding equilibrium result calculated from 
using £o = V ( \2^ ). In addition to the overall smoothing 
of the singularities, expected in a non-equilibrium system, 
there are two significant non-equilibrium features. Firstly, 
the maximum in the spectral weight is shifted away from its 
equilibrium value by an amount proportional to the applied 
voltage. The shift, Re(A(V) — A(0)), which is given in the 
forward- biased case in 1B4L depends on the polarity of the 
voltage. Secondly the form of the function changes on revers- 
ing the polarity of the device. 



—i times the logarithms of the eigenvalues of the scat- 
tering matrix S e ), which appear in the formulas (|53[l . 
are replaced by complex phase shifts. In the forward 
bias case described by l|53ll . these are — ilogS*^ and 
— zlog(l/S , | 2 )*- The effect of the complex phase shifts is 
to smooth the singularity seen in equilibrium (this could 
be expected on quite general grounds) and to introduce 
a polarity dependence. This polarity dependence affects 
both the shape and the position of the spectrum and is 
evident in Figure [3] where we show xd^) for a particular 
choice of S e . The dependence of the spectrum, xg(u>), on 
the polarity of the device, when operating out of equilib- 
rium, is governed by the difference «i2 = ot\ — a-i (with 
ai.2 as defined in the figure caption). The difference in 
the overall position of the spectrum on changing the po- 
larity is given by the difference in the second term on 
the right hand side of (|65() and is proportional to ot\2- 
This origin of this shift of the spectrum is the change in 
the nature of the scattering across the barrier from fully 
coherent in the equilibrium case to incoherent for times 
tf ^> V^ 1 in the non-equilibrium case. The shape of 
the spectrum reflects the decay of charge from its initial 
distribution (the equilibrium distribution for S = 1) to 
the steady-state distribution for S = S e 0, . In the 
non-equilibrium case, this decay can occur differently de- 
pending on the polarity. If more charge is needed in the 
left-hand electrode to screen the potential characterized 
by S e , than in the right-hand one (a±2 > 0), this charge 



FIG. 4: Spectral function KexF{^) computed from IBlOl with 
ui in units of the bias voltage V for the case (Wi, W%) ~ (0, 1) 
in 1511 1. The curve marked + (— ) refers to the case in which 
electrode 1 (2) is at the higher chemical potential. Also shown 
is the corresponding equilibrium result [gsj. For these rela- 
tively small phase shifts the singularity seen in equilibrium 
disappears completely, although there is still a polarity de- 
pendence of the spectrum even though the scattering matrix 
is symmetric. 



can come from states within V of the Fermi energy of 
the right-hand electrode when the device is reverse-biased 
but not when it is forward-biased. 

For the model device shown in Fig [3 the absorption 
spectrum is given by the Fourier transform (see Appendix 
FBI of G(tf) rather than F(tf), as the transition in the 
barrier is presumed to be dipolar. However, the corre- 
sponding function F(tf) is also important. In HcHjHJ Yu- 
val and Anderson showed that the Kondo problem could 
be treated as an infinite sequence of spin flips or switching 
events, with the response of the conduction electrons to 
each switching event characterized by F(tf), Given the 
long-standing interest in non-equilibrium effects in the 
Kondo effect [32, , the correct non-equilibrium form 
for F(tf) would be the starting point for the study of 
the non-equilibrium Kondo effect using a generalization 
of the Yuval- Anderson mapping. 

We must first write the function F and the open line 
function L in terms of the gauge-transformed basis: 



L(t f ) ~ W*^e- iLVt t 



(68) 

We can now insert the solution for Y + (with A = 1) from 
(|61fl into 151fl . The result can be written: 



L{t f ) 



it f 



-W 



[az f (iVt f ) 2 *+ 



W, (69) 



where x± = xi(l) ± £2(1) with xdX) defined in 16311 . 
a = S^/Sf-t and z/ = e - iVt f 34]. The absence of a 
contribution proportional to W^W2Zf is to be expected. 
This would involve a contribution to the open-line func- 
tion from an electron initially placed in the right-hand 
electrode exciting the Fermi gas in the left electrode. 



11 



Since we are assuming that the tunneling through the 
barrier by the electron is a slow process on the scale 
of 1/V, this does not lead to a singular contribution to 
F. (There is still a contribution to F proportional to 
W*W2Zf from the direct scattering term CG(tf) in 1481 ) 
The effect of the open line contribution on F{ti) = 
G(tf)L(tf) is the natural generalization of the equilib- 
rium result that one might expect given the results for 
G(tf). The corresponding spectral functions is shown in 
Figure 0] for a particular choice of S e . For simplicity we 
only look at the case where the electron is added and 
removed from the same (i = 2) electrode, ie W ~ (0, 1). 
The dependence of the spectrum on the polarity of the 
device, is present even in the case where S e is symmetric. 
When a particle is added to an electrode, the response of 
the system will be depend on whether the electrode is at 
the higher or lower chemical potential. The form of the 
spectrum can also differ substantially from what happens 
in equilibrium. For the case a.\ = ai = 0.4 and R = 0.7 
shown in Figure^] there is no real peak left over from the 
equilibrium result. This is because the phase shifts 5^, 
corresponding to the eigen modes of S e , are small, and 
hence the exponents in (|53[1 for the equilibrium function 
Xf(w) ~ u}~( s <~ v > /*■ are also small. The corresponding 
singularity is weak and easily smoothed out by the finite 
lifetimes of states close to the Fermi energy in the non- 
equilibrium case. This smoothing is enhanced because 
one of the phase shifts, 5^, is always larger than ot\ and 
Q'2- This larger exponent gives the dominant singularity 
in equilibrium, but is then effectively replaced by a\ out 
of equilibrium. 



IV. CHARGE TRANSFER: (CSAC) 

The existence of Coherent States of Alternating Cur- 
rent (CSAC) was predicted in [[J. These consist of a 
sequence of pulses which propagate through a contact. 
When the bias across the contact is described by a class 
of periodic (with period fi) rational functions of the vari- 
able z = e lSlt , then the shot noise is minimised and the 
noise distribution does not depend on the separation of 
the pulses. This result is still not well understood, nor 
is it possible, using the original derivation, to establish 
how robust these states are against deviations from zero 
temperature or from the ideal pulse shape. 

Recent rapid experimental progress in the application 
of micro-wave radiation at low temperatures suggests 
that the experimental test of the existence of the CSAC 
is just about possible. Several experimental groups are 
pushing the technology in this direction [3^, |3(| |3g] , and 
it should only be a matter of time before experimental 
data becomes available. However, interpretation of these 
future experiments will not be easy using the analyti- 
cal method used in 0] as this depends crucially on the 
particular shape of the pulses. There are no predictions 
about what happens when the shape of the pulses devi- 
ates slightly from the required one (something unavoid- 



able in any real experiment), nor is the effect of non-zero 
temperature known. 

Here we show that the results of pj for the CSAC are 
easy to derive using the RH approach. When the pulses 
are periodic as in the case of the CSAC, the RH problem 
simplifies. It requires solving for functions which are an- 
alytic in two disconnected regions (inside and outside 
the unit circle) with the jump function specifying the dis- 
continuity across the boundary between them. We show 
that the particular case of the CSAC corresponds to an 
RH problem which can be solved exactly using combi- 
nations of meromorphic functions in the plane — one of 
which is analytic inside and one outside the unit circle. 

The model device considered in 0,0 consists of a tun- 
nel junction driven by a bias voltage V(t) which is pe- 
riodic in time with period T = 2tt/Q. It is equivalent 
to the device shown in Fig [2 We are interested in the 
change in physical quantities over one cycle of the pump 
in the limit tf 3> T. In this limit, effects induced by the 
switching on and off of the periodic potential at t = 
and t = tf are irrelevant. The scattering matrix S is at 
some constant value S e between and tf. Applying the 
time-dependent gauge transformation leads to S be- 
coming a periodic function of time, so that it no longer 
commutes with itself at different times. 

The distribution function for any single-particle ob- 
servable measured in this periodically pumped Fermi sys- 
tem will involve the solution of a non-commuting RH 
problem. In particular, the characteristic function or gen- 
erating function for moments of the distribution of the 
net transfer of charge from electrode 1 (left electrode) 
to the electrode 2 (right electrode), x(A), is given by (0} 
with the operator R given by 

R(\) = U^{tf)er lX ^^il{tf)e lX ^. (70) 

Here 

Q^Y.^La,. (71) 

c 

For states close to the Fermi surface (E = e + e' ~ 0), 
the matrix R in the time-representation can be written 

R(t, A) = S-\t)e- iXL S(t)e iXL , (72) 

so that the characteristic function will be given by (|12l 

1131113 

Iog X (A)=Tr(ln(l-/ + /i?)). (73) 

If the inverse of the solution Y + (t) to the RH problem 
(|22l I23fl with R given by l|72|) , is analytic in the upper 
half-plane, we can write the characteristic function as 

logx(A) = J dX' J dttt^Y+fY^R- 1 ^} (74) 
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Using H72f> and (|21|l , and computing explicitly the deriva- 
tive with respect to A, we obtain |2ll38|| 



rA d\' 



dttr ■ 



(It 



{Se lA L Y + 



(75) 

(If the eigenvalues of Y+ have zeros in the upper half- 
plane, there are additional contributions to the right 
hand side of the corresponding relations to (|24|) for Y7 1 
from its poles. In this case, (1 — / + fR)^ 1 is no longer 
given by (|25|) but can be found using methods described 
in Chapter 6 of 

In the case of the periodically driven pump, the scat- 
tering matrix (after applying the gauge transformation) 
is periodic S(t) = S(t + T). If we change variables to 
, we need to find a function Y" + (z), which is 



a 27Tit/T 



analytic for \z\ < 1 and Y"_ which is analytic for \z\ 
and Y"_ — > 1 when \z\ — ► oo. On the unit circle \z\ = 



Y-Y+ 1 = S- 1 (z)e- tXL S(z)e 



iXL 



> 1 

= 1, 

(76) 



The characteristic function for charge transmitted during 
one cycle of the periodic pump in the limit tf/T^> 1 is 
given by 



log'X 



dz 



tr 



dX 

^ J\z\ = l 

d (Se lXL Y + 
dz 



(77) 



i\L 



Y- 



L 



CSAC's were reported in 0] for the case when the 
phase factor in l|58(l can be written as a rational func- 
tion, l(z), of the variable z — e 27Tlt / T : 



,+iL J* V(r)dr 



l(z) 

N 

'(*) = Ut 



i=l 



ajz 



(78) 



where either all |oi| > or \cii\ < 0. We can choose \<n\ > 
without loss of generality as z i— * l/z simply reverses 
the polarity of the device. In this case, we decompose 
R(t, A) (see @) as follows 



R = 



1 

q 



a 
i 



1 Pl(z) 
1 



(79) 



where 



|S 12 |V A + |S n | 

^21 (^22) (l " 







(80) 



The solution to the RH problem 
1 

a 



-1 



a 
i 



1 

1 



is clearly 



1 
1 



7R 




1 



a 



Inserting this into l(77Jl gives the result reported in 0: 



ln X (A) 



N^hxa 
2tt 



(82) 



with a given by l|80() . 

The surprising feature of the result l)82|l is that it im- 
plies that the second moment of the shot noise <n 2 > 
achieves the absolute minimum for given charge transfer 
< n > 0, which is the value obtained in the constant 
bias case — > 00 for all i. The feature of the phase factor 
l(z) which leads to the RH problem being so easy to solve 
is that all its poles (zeros) are either inside or outside the 
unit circle \z\ = 1, which means that the decomposition 
of R in (|79(l automatically solves the RH problem. In the 
case of an arbitrary rational function for l(z) this is not 
the case as there can be points at which dct|F + (z)| van- 
ishes inside the unit circle. The corresponding formulae 
to (|24fl Y+ 1 pick up additional terms on the right hand 
side nad (1 — f + .fR)~ 1 is not given by l|25|). although, in 
principle, it can still be found given the solution to the 
RH problem Y(z). 



V. CONCLUSIONS AND OUTLOOK 

The RH approach is a general method for computing 
the response of a Fermi gas to a localized time-dependent 
perturbation. There are two key steps to the method. 
First, provided the condition J5J is met, the method 
works with the scattering matrix defined on the instanta- 
neous value of the potential rather than with the Hamil- 
tonian. This has the attractive feature of working di- 
rectly with the physical quantities determining the long 
time response of the system to a perturbation, namely 
scattering amplitudes for particles close to the Fermi sur- 
face. The condition J5J is essentially the requirement that 
the perturbation varies more slowly than the delay time 
for a particle traversing the region in which the pertur- 
bation acts. The second key step is to relate the response 
of the Fermi gas to the solution of a non-commuting RH 
problem 1221 123|) . The RH problem corresponding to any 
given experimental situation is usually easy to set up. Its 
solution and the interpretation of the results is a more 
delicate task which needs to be repeated for each new 
physical situation. While there is no analytical solution 
of the general non-Abelian RH problem, there is a pow- 
erful technique for finding asymptotic solutions valid for 
frequencies much smaller than those present in the jump 
function [j]. 

Here we have emphasized the generality of the ap- 
proach and applied it to two existing problems — the 



13 



Fermi Edge Singularity and the shot noise in a period- 
ically pumped tunnel junction. The calculations in the 
two cases are very similar. In the case of the FES we have 
rederived all the known results for the equilibrium case 
emphasizing, in particular, how the method is no more 
complicated in the case of the non-separable potential 
than in the separable case. Our derivation for the non- 
separable case is, we believe, the first non-perturbativc 
solution of this problem. For the non-equilibrium device 
shown in Fig|2 we have explained how the results for the 
core-hole Green's function of were obtained and given 
the corresponding results for the open-line function L(tf) 
(|69|1 . For the case of the CSAC's, we have shown that 
the particular form of the periodically varying bias with 
the phase factor Z(e 27 "*/ T ) given by 1)78(1 corresponds to 
a case in which the RH problem can be solved exactly. 

It is possible within the RH approach to handle correc- 
tions to the asymptotic solution to the non-communting 
RH problem we have been using in order to allow us 
to compute the response of systems in the intermedi- 
ate regime (where one is interested in the response at 
frequencies comparable to those introduced by the per- 
turbation). The RH problem lends itself naturally to a 
type of perturbativc analysis. The corrections to the ap- 
proximate solution valid for long times, <|f>()(l and 
can be described by multiplying the approximate solu- 
tion by a function which is analytic except across the 
additional vertical cuts introduced to simplify the origi- 
nal problem. This function can be specified by a Cauchy 
integral around the cut. Preliminary work in this direc- 
tion has been attempted in [^j . 

Finally, the RH method should generalize to non-zero 
temperatures. As was observed in [2(|, the singular in- 
tegral equation appearing at finite temperatures in a re- 
lated problem can be solved analytically. Also, the ana- 
lytic treatment of the finite-temperature Fermi-edge sin- 
gularity in 0, |4£j again suggests that the RH approach 
will generalize successfully to finite temperatures. 



APPENDIX A: BOUND STATES 

If the perturbing potential generates a bound state(s), 
then © is no longer correct. In the case where the po- 
tential (and hence S(t)) simply switches between its un- 
perturbed value and a new but time-independent value 
at t = and back again at t = tf, we can correct a by 
including the effect of the bound state explicitly. The 
treatment follows closely that of ^jj> although only the 
case of a separable potential was treated there. We write 



+ e iHot \b}e- iEbt (b\. 



(Al) 



Here \b) is the bound state wavefunction, while a de- 
scribes the scattering of the states within the continuum, 
and is given by the Fourier transform of the scattering 
matrix S(t) @) as before. (Hq is the matrix of Hq taken 
between single-particle basis states.) 



For the case of the function G(tf) = det|l — / + fa\ 
(see I31fl we have 



G(t f ) = G(t f )det\l + A\b}(b\\ 
= G(t f )(l + A B ), 

where A B = (b\A\b) with 

A = (1 - / + fa)- l fe iHot e- iBbt , 



(A2) 



(A3) 



and where G(tf) = det|l — /+ /cr|. We write the bound 
state as an expansion over the basis vectors 



(A4) 



For long times tf the response is dominated by states 
within 1/tf of the Fermi energy and it is a reasonable ap- 
proximation to neglect the energy dependence of the co- 
efficients u c . After switching to the time-representation, 
and using (|49|l with a in place of a, we obtain [yi is the 
density of states in channel I) 

(b\A\b) = e - iEbt f [ dede ' [ dhdt-iUpUi x (A5) 

/ ie't! 



Integrating over energies and times gives 



A B ~ e- m ^ Ul JV x 



(A6) 

Provided S(t) commutes with itself at all times between 
and tf, Y is given by (|37() . For the single-channel case 
with S = e 2lS , we obtain 



(b\A\b) 



(A7) 



Here we introduce the quantity S which is the phase shift 
modulo 7r and takes values on the interval [— ir/2, 7r/2]. 
Normally the phase shift S is defined with a jump of it 
at a bound state thereby ensuring compliance with the 
Friedel sum rule pol |. However, when writing the scatter- 
ing matrix as in IjAlfl . the contribution from the bound 
state to the scattering matrix is explicitly included in the 
second term on the right hand side and is not in the scat- 
tering matrix S. At the bottom of the band, the value of 
the phase shift which enters the threshold shift is clearly 
5 as emphasized in |llj . 

Although the calculation is longer, the function F(tf) 
can be obtained in a similar manner by replacing a in 
(|46|l by the form lAlfl . One needs only to keep track of 
terms up to first order in e ~ tE b t f _ (Higher order terms 
must give zero as they correspond to double or higher 
occupancy of the bound state. They can be seen to make 
no contribution by subsituting the formula <A1|) in 1421) .') 
As for the case of the function G(tf) considered above, 
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we neglect the energy dependence of W e and u £ (see l2"51 
and lA4|l . We define 



C b = (W* ■ u)(u* ■ W)e 

o = (l-z + j^r 1 /. 



-iE b t f 



(A8) 
(A9) 



Here O is just the scattering state contribution to O (see 
123: 

O ~ 6 - Oe lHat * \b)e- iEht f (b\0. (A10) 

We obtain (from EH} 

F(t f ) = G(t f )[C-(g\0\h)\ 

~ G(t f )C b - G(t f ) [1 + A B ] (g\0\h). (All) 

Retaining the dominant terms and ignoring the possi- 
bility that there is an unexpected cancellation between 
terms proportional to e _i£6 *^ , 

F(tf) = G(tf)L(t f ) + aC b G(t f ) (A12) 

where a ~ 1 is some constant and L(tf) is the scat- 
tering state contribution to the open-line function. For 
the single channel case with S = e l2S , we again assume 
that the exponent in G(tf) is 5 = 5 — n and obtain 
F(t f ) ~ F b (t f ) + F {tf) with 



F b (tf) 



-iE b t, 



(iS)*/)* 4 /*- 1 ) 3 
(A13) 



APPENDIX B: COMPUTING SPECTRAL 
FUNCTIONS 

Given G(tt) or F(tf) we would like to compute the 
corresponding spectral functions given by a Fourier in- 
tegral over tf. Assume that scattering matrix, S e , has 
diagonal elements VRe i2ai - 2 . Using the complex cutoff 
(V (normally £ = i), we have from (|64|) 

logG(f /s V) = -i(E - A(V))t f - p G log (CVt f ) + D. 

(Bl) 

The exponent Pq = x\ + x 2 _ , where 

x+= \ogSf 1 ^a 1 _AogR 
2iri 7r 47r 

and 

g _ = log (1/^5) = o i+ .lggfl j 

27TZ 7T 47T 

The modified threshold shift is given by (|65|l 

A(V) = f A(0) + V ^^S%A _ lV (]^ 

\ 7T / \ 47T 

(B4) 



The real part of A(l/) fixes the threshold. We will ab- 
sorb this into the definition of frequency when computing 
Fourier transforms. 

We write 0c = Pgi + i@G2 with 



0G1 = 



«2\ 2 



logi? 



2jt 



and 



CI 



G2 



(ai - a 2 ) logR 
7r 2ir 



(B5) 



(B6) 



For the function F(tf), the exponent becomes (3f — 
(x- - l) 2 + x\ or j3 F = x 2 _ + (x + — l) 2 depending 
on whether the electron is added to the electrode with 
lower or higher chemical potential. This gives f3p = 
Pa — 2x± + 1 and 

Pfi=Pgi-^ + 1 and f3 F2 = f3 G2 ± ^ . (B7) 



2vr 



Introducing 



UJ 2 = -l0gi?/47T, 



(B8) 



the spectral functions of G or F are proportional to the 
real part of the Fourier integral, Xf,g{c)i where: 



X(e) 



dt f (cvt f ) 



-0 (ie-uj2V)tf 

i 



(B9) 



with (3 = (3f for Xf and Pg for XG- Here the lower limit 
of the integral is taken to be 0, which is only valid when 
Pi < 1. When Pi > 1, contributions from the lower limit 
of the integral dominate and the response is dominated 
by high frequency contributions which are not changed 
from the equilibrium case. These are not described by 
the formula l|64(l and depend on details relating to the 
band edge. If the phase shifts x± are small, which can 
be the case for the spectral function of F (or GAb in the 
presence of a bound state, see IA2I and IA7|) . then Pi will 
be close to 1. In this case x(e) given in l)B9(l contains a 
significant contribution from times tf < 1/V for which 
our asymptotic solution for F(tf) (or GAb) is incorrect. 
We can correct for this by noting that when 1 — /3 <C 1 
the contribution from times with Vtf < 1 gives just a 
constant offset which can be subtracted from \- To see 
this, we expand the exponential term e ( ie -"2V)t j n ^ e 
integrand and integrate term by term from Vtf = to 
Vtf = 1 . The first term in the expansion is independent 
of e and much larger than subsequent terms provided 
(e/V) <C 1/| 1 — P\- In practice we subtract from the real 
part of x its value at u> ~ —V. (When 1 — P is not small 
the contribution from the times Vtf < 1 to the real part 
of x is negligible anyway. ) 

Eq l|B9|l is in the form of a standard integral and (see 
8.312.2 in [HI) is given by: 



i 



V \ wi + iuj 2 



1-0 



T(l-p). (BIO) 
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If we define Q = \uji + ici^le 1 ^" and write: The real part of xi^i) carL then be written p9l |: 

r(l-/?) = |r(l-/3)|e^ and iC = e*+<, (Bll) \V(1 - 8)\ 1 

then v n Pl 

sin [/3i(0 c - 4>n) + (<fo - fa) - fh logfi]. 

• „i/3 2 log n 

x(o>i) = e - iM < +M < — - x 

V Q 1 ^^ 1 For both functions F and G, the cutoff parameter ( = i, 



J(Pi-i)^-Mo\ r ^ _/3)| e ^r. (B12) 



SO 0/- = 7T. 
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